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Abstract 

MgGeOa-perovskite is known to be a low-pressure analog of MgSiOs-perovskite in 
many respects, but especially in regard to the post-perovskite transition. As such, in¬ 
vestigation of spin state changes in Fe-bearing MgGeOs might help to clarify some 
aspects of this type of state change in Fe-bearing MgSi 03 . Using DFT-hU calcula¬ 
tions, we have investigated pressure induced spin state changes in Fe^”^ and Fe^^ in 
MgGe 03 perovskite and post-perovskite. Owing to the relatively larger atomic size of 
germanium compared to silicon, germanate phases have larger unit cell volume and 
inter-atomic distances than equivalent silicate phases at same pressures. As a result, all 
pressure induced state changes in iron occur at higher pressures in germanate phases 
than in the silicate ones, be it a spin state change or position change of (ferrous) iron in 
the perovskite cage. We showed that iron state transitions occur at particular average 
Fe-O bond-length irrespective of mineral composition (silicate or germanate) or func¬ 
tionals (LDA-i-Ujc or GGA-hUjc)- Ferrous iron substitution decreases the perovskite 
to post-perovskite (PPv) transition pressure while coupled ferric iron substitution in¬ 
creases it noticeably. 
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1. Introduction 


Iron-bearing MgSiOs perovskite (Si-Pv) is one of the major constituent minerals 
of Earth’s lower mantle along with (Mg,Fe)0 ferropericlase (Fp). Unraveling the com¬ 
position and thermal structure of the lower mantle requires a detailed understanding of 
the influence of iron (Fe) on these host minerals. It is well known that iron in (Mg,Fe)0 
undergoes a pressure induced spin crossover from high-spin (S=2) to low-spin (S=0) 
state and affects elastic properties ^adro et al.\ |2003| |2004[ \Goncharov et al.\ |2006[ 


Cwwhurst et al. 2008 [ Tsuchiya et al. 2006 [ Marquardt et al. 2009[ Wu et al. 2009 


\Antonangeli et g/TlpOl 1 [[Wm et Q/.]|2013[|7f5M and Wentzcovitch\\20l^\Wu and Wentz-\ 

covitch\ 2014|l. In the case of Si-Pv, deciphering the iron spin state has been quite 


challenging due to complex perovskite structure and possibility of different valence 
states of iron (ferric and ferrous). Si-Pv is also known to undergo a perovskite (Pv) 
to post-perovskite (PPv) phase transition at ~ 125 GPa and this will have geophysi¬ 
cal consequences in deep lower mantle region ^Murakami et al. \ 2004[ Tsuchiya et al.\ 
2004[ Oganov and Ono\ 2004|l. Understanding the effect of iron, especially, on Pv to 


PPv transition and on elastic properties of Pv and PPv phases is crucial to constrain 
lower mantle composition. 

In (Mg,Fe^^)Si03-Pv, the pressure induced increase in Mdssbauer quadrupole split¬ 
ting (QS) above ~30 GPa (low-QS to high-QS transition) [Jackson et al. 2005 [ Li etal.\ 
2004| 2006[ McCammon et al. 2008 |2010[ Lin et al. 2012 |2013| l has been explained 


by lateral displacement of Fe^^ ion which remains in the high-spin (HS) state in the 
entire lower-mantle pressure range ^Bengtson et al. 2009[ Hsu et al. 2010[ Hsu and 


Wentzcovitch 20141. In the case of (Mg,Fe^^)(Si,Fe^^)03-Pv, however, hrst-principles 

(201 1| repotted that Fe^^ in B-site undergoes a high-spin 


calculations by 


Hsu et al. 


(S=5/2) to low-spin (S=l/2) state change at approximately 41 GPa and 70 GPa using 
LDA-hUjc and GGA-hU^c methods, respectively, while in the A-site Fe^^ remains in 
the high-spin state. These findings were in good agreement with experimental obser¬ 
vations of \Catalli et al. ( |2010a I. There have also been several experimental ^Lin et aL\ 
2008[ Catalli et al. 2010b Mao et al. 2010|l and first-principles [Zhang and Ogano^ 


2006[ Caracas and Cohen 2008 Yu et al. 20121 studies about the spin crossover in 
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MgSiOs-PPv (Si-PPv). The consensus regarding the spin crossover in iron-bearing Si- 
PPv is that Fe^^ in the A-site remains in the high-spin and high-QS state throughout 
the lower mantle pressure range, while Fe^^ undergoes a high to low spin transition in 
the B site and remains in the high spin state in the A-site \Yu et al. 2012| l. In spite of 
these efforts to pin down the state of iron at high pressures, it has been very difficult 
to ascertain the exact role of iron with its varying spin and valence state on Pv to PPv 
transition, especially because of the extreme pressure and temperature conditions in the 
deep lower mantle region. 


MgGeOs-perovskite (Ge-Pv) is a low pressure analog of MgSiOs perovskite {Ross 


and Navrotsky 1988| Leinenweber et al. 2004[ Hirose et al. \ 2005 1 Kubo et al. 2006| 

\Merkel et al.\ |2006[ \Runge et al.\ |2006[ \Shim et al] |2007[ \Tsuchiya and Tsuchiya\ 
2007[ Duffy\ 2008 Ito et al. 2010|l and would be a promising candidate to investigate 


iron spin crossover and its structural influence on MgSi 03 -Pv and PPv phases. In this 
study, we have used first-principles DFT-i-U (LDAh-Ujc and GGA-hU^c) calculations to 
investigate the iron spin states and their crossover in the Pv and PPv phases of Fe^^- 
and Fe^^-bearing MgGe 03 . Comparing enthalpy of the different spin states, we have 
determined the stability of the different states and transition pressure. Using these 
results, we have also investigated the effect of Fe^’^ and Fe^^ substitution on Pv to PPv 
transition. 


2. Method 

We have used density functional theory (DFT) within the local density approxima¬ 
tion (LDA) ( Ceperley and Alder 1 19801 1 ^nd generalized gradient approximation (GGA) 
{Perdew et al. 1996| l augmented by the Hubbard U (LDA/GGA-t-Ujc), calculated self- 
consistently and structurally consistently ( Cococcioni and de Gironcoli 2005 [ Kulik 
et al.\ 20061 2009 [ 2011} Himmetoglu et al. 2014| l. Ultrasoft pseudo¬ 

potentials, generated by Vanderbilt’s method ^VanderbJii\ |1990| l, have been used for 
Fe, Si, Ge, and O. For Mg, norm-conserving pseudo-potential generated by von Barth- 
Car’s method, has been used. The plane-wave kinetic energy and charge density cutoff 
are 50 Ry and 300 Ry, respectively. Self-consistent Ujc, calculated by Hsu et al. ( |2010 1 
and Yu et al. (2012|l using the linear response approach ( Cococcioni and de Gironcoli 


3 































































































2005 [ Kulik et al. 2006| l (Table[T]l, have been used. All calculations were performed in 
a 40-atoms super-cell. Electronic states were sampled on a shifted 2x2x2 and 3x2x3 
k-point grid for Pv and PPv structure, respectively ^Monkhorst and Pa^ 1976[ ). Struc¬ 
tural optimization at any arbitrary volume has been performed using variable cell-shape 
damped molecular dynamics approach ^Wentzcovitch l99l\\Wentzcovitch et al. 1993|l. 


Structures are optimized until the inter-atomic forces are smaller than 10 Ryfa.u. For 
every possible spin and valence state, structures are optimized at twelve volumes in the 
relevant pressure-range and the 3'''^-order Birch-Mumgharn equation of state has been 
used to fit these compression data. 


3. Iron spin states in Pv and PPv phases 


Iron spin states have been investigated in Fe^^ and Fe^’^-bearing MgGeOa and 
MgSiOs in a 40-atom super-cell. One Fe^^ ion is substituted for a Mg^^ ion in the A- 
site in (Mgo. 875 Feo.i 25 )Ge 03 and (Mgo. 875 Feo.i 25 )Si 03 while two Fe^’^’s are substituted 
for one Mg^^ and one first neighbor Ge"^^ or Si"^^ as in (Mgo. 875 Feo.i 25 )(Geo. 875 Feo.i 25)03 
and (Mgo.875Feo.i25)(Sio.875Feo.i25)03. The nearest-neighbor Fe^^-Fe^^ substitution has 
been chosen as it has been found to be the lowest energy configuration in previous first- 
principles studies ^Hsu et aT\ |201 1| |2012[ \Yu et al.\ |2012| l. Relative enthalpies have 
been used to calculate transition pressures (P,) and the stability of a state or phase in a 
given pressure range. Effects of disorder have not been addressed at this point, but it 
is expected to produce a two-phase loop. Here we wish to analyze (a) the potentially 
different effects of Fe^^ and Fe^^ in the Pv to PPv transition and (b) whether the ger- 
manate phase could serve as a good low pressure analog to investigate spin transitions 
in iron. 


3.1. Spin states of in Ge-Pv 

As in the case of Fe^’^in Si-Pv [Bengtson et 
Wentzcovitch 2014 1 , Fe^’^ in Ge-Pv remains in the high-spin (HS, S=2) state in the 


2009 

Hsu et ail 

2010 

Hsu and 


entire lower-mantle pressure range (0-135 GPa). However, it undergoes a pressure in¬ 
duced lateral displacement producing a transition from low to high Mdssbauer quadrupole 
splitting (QS). The atomic structures for the low- and high-QS states of iron in Si-Pv 
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are shown in Fig. [Ja) indicating the lateral displacement of iron in the perovskite 
A-site. 

The pressure dependence of enthalpy differences between high-QS and low-QS 
states {AVi-Hhigh-Qs-i^iow-Qs) gives the low- to high-QS transition pressure. It is well- 
known that LDA usually underestimates transition pressures while GGA usually over¬ 
estimates it. In this study, we present both LDA-rUic and GGAh-Ujc results to provide 
lower and upper bounds for this transition, as shown in Fig. [TJb). The low- to high- 
QS state change in Ge-Pv occurs at ~28.5 GPa with LDA-hUjc and at ~40.5 GPa with 
GGA-rUic. These transition pressures are quite large compared to those in Si-Pv, i.e., 
~9.5 GPa and ~22.5 GPa respectively as presented in Table]^ 

Owing to the larger ionic radius of Ge^^ (0.53 A) compared to that of Si"^^ (0.40 
A), the unit-cell volume and inter-atomic distances in Ge-Pv are larger than in Si- 
Pv. Zero pressure Fe-O bond-lengths in Si-Pv and Ge-Pv obtained with LDAh-Ujc 
are shown Fig. |^a). The average of Fe-O bond-lengths, (Fe-O), is also shown in 
Fig. I^b). The larger (Fe-O) value in Ge-Pv allows Fe^"^ to sustain higher pressures 
than in Si-Pv before changing to the high-QS state. It is also interesting to note that 
this transition occurs when the (Fe-O) reaches a particular value (~2.22 A) whether 
in Ge-Pv or in Si-Pv, whether using LDA-rUic or GGA-hUjc methods (Fig. |^(b)). It 
should be noted that the Pv to PPv transition is expected to occur at approximately ~63 
GPa in MgGe 03 \Hirose et a/.||2005] l and at ~ 125 GPa in MgSi 03 ^urakami et al.\ 
2004^ . These pressures are much higher than low- to high-QS transition pressures in 
both Fe^^-bearing compounds. Consequently, this transition has not been investigated 
in Si-PPv or Ge-PPv phases. 

3.2. Spin states of in Ge-Pv and -PPv 

Next, we investigate spin states of Fe^^ in the B-site of Ge-Pv (low-pressure) and 
PPv (high-pressure) phases. The crystal structures of Si-Pv and Si-PPv phases are 
shown in Fig.j^a). Si-PPv has a layered structure characterized by corner-sharing and 
edge-sharing SiOe octahedra extending along the a-c plane. The pressure dependence 
of calculated enthalpy differences between high-spin (HS) and low-spin (LS) states 
of Fe^^ (AH=I-1 b(/,s)-I-1b(hs)) in Pv and PPv phases are shown in Fig. [^(b) and (c). 
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respectively. The calculated HS to LS transition pressures in Ge-Pv are 56.5 GPa with 
LDA+Ujc and 85.0 GPa with GGA+Ujc. Similarly, these transition pressures in Ge- 
PPv are 49.5 GPa and 82.5 GPa with these methods. Our calculated transition pressures 
in Si-Pv are 41.0 GPa with LDA+Ujc and 69.5 GPa with GGA+Ujc, which are in very 
good agreement with previous first-principles study ^su et oT. 201 1| |. In the case of 
Si-PPv, our calculated transition pressures are 28.0 GPa and 59.5 GPa, respectively, 
which are quite larger than previously reported values of -30 GPa and 36 GPa ^¥u\ 
et al.\ 2012|l. These transition pressures in Si-Pv and Si-PPv are also very similar to 


those reported by Hsu et al. (20121 in (Mg,Al)(Si,Fe) 03 -Pv and PPv phases where it 
was observed that the presence of aluminum does not affect the spin transition pressure 
significantly. Considering the overall agreement of the stability of the iron spin states 


and pressure induced crossover in Si-Pv and Si-PPv phases with previous studies (Hsu 


et al. 2010|| 2011 20121, the present calculations about HS to LS transition in Fe^^ at 
B site seems to be consistent and reliable. 

In order to explain the increase of HS to LS transition pressures in MgGeOs phases, 
we compare atomic configurations around Fe^^ in the B-site of Si-Pv and Ge-Pv phases. 
As shown in Fig. |^a), the coordination number of Fe^^ is six and Fe-O bond lengths 
in this site are smaller than those of Fe^^ in the A-site. The pressure dependence of 
the average bond-length, (Fe-O), for Fe^^ in the B-site in Si-Pv and Ge-Pv phases are 
shown in Fig. 0b). As in the case of low- to high-QS transition in Fe^^ (Fig. [^, the 
HS to LS transition in Fe^^ in Ge-Pv and Ge-PPv phases are also large compared to 
those in Si-Pv and Si-PPv phases (Fig. due to relatively larger Fe-O bond-lengths 
and unit-cell volumes in germanate phases at a given pressure (Fig. |^. Similar to low- 
to high-QS transition, the HS to LS transition also occurs only when the (Fe-O) reaches 
a particular value (i.e., ~L86 A in the present case) irrespective of Ge-Pv or Si-Pv, or 
of the exchange-correlation functional used (Fig. I^b)). This is a useful insight for 
predicting pressure induced spin state changes in different materials. 

4. Effect of Fe^"^ and in the Pv to PPv transition 

Having investigated the stability of different iron states in the lower mantle pressure 
range in the previous section, we now focus on the potential effect of Fe^’^ and Fe^’^ 
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substitution on the Pv to PPv transition. This is a plain static calculation to access the 
overall influence of iron substitution on the Pv to PPv transition in these phases. We 
have not included finite temperature entropic or vibrational effects, therefore, we are 
not able to produce a full phase diagram with a two phase loop region. This is simply 
the first step towards the calculation of such phase diagrams. 

It is well known that the LDA functional underestimates the equilibrium volume 
while GGA overestimates it. After inclusion of vibrational effects LDA equations of 
state and elastic properties tend to agree better with experimental data than GGA results 
\Karki et al.\ |1999t \Karkiet gT] |2000[ \Wentzcovitch et aT\ |2004[ \Nunez-Valdez et al] 

2012a|b 2013| l. Therefore, we have chosen LDA+U^c results to investigate the effect of 
Fe^^ and Fe^^ in the Pv to PPv transition in both compounds. Previous studies based 
on ab initio calculations \Tsuchiya et al. \ |2004[ \Tsuchiya and Tsuchiya\ |2007| l and 
experimental measurements ^Murakami et Q/.[[2004[|//i>o5e et q/.|| 2005| | have reported 
the Pv to PPv transition in MgSi 03 and MgGeOa to occur around ~125 GPa and ~63 
GPa, respectively, at ~2500 K. 

As indicated in Fig. [TJb) and in Table the low- to high-QS transition in Si- 
Pv and Ge-Pv occurs at much lower pressures than the Pv to PPv transition in these 
minerals {Murakami et al. 2004 Hirose et al. 2005 | l. Considering this fact, only the 
high pressure high-QS state of Fe^^ was considered in the investigation of the Pv to 
PPv transition. Similarly, near the Pv to PPv transition pressure, the stable state of 
Fe^^ substituted in both Si-Pv and Si-PPv phases is HS state in the A-site (HS-A) 
and LS state in the B-site (LS-B). However, in germanate phases, the HS-B to LS-B 
transition pressure (Fig. and Table is near the Pv to PPv transition pressure. In 
order to find the appropriate spin states of Fe^^ across the transition we have considered 
all possibilities in the calculation of enthalpy differences (AH=Hpp,,-Hp,.) (see Fig. |^. 
As indicated in Fig. |^b) the transition pressure 55.5 GPa involves a transition between 
HS-B in Ge-Pv to LS-B in Ge-PPv. Obviously iron in the A site remains in the HS state. 
At high temperatures, contributions by other states might be relevant as well, but this 
issue is well beyond the scope of this paper. This transition and transition pressure, i.e., 
HS-B in Pv and LS-B in PPv, is the only combination consistent with B(HS) to B(LS) 
transition pressures of 56.5 GPa and 49.5 GPa in Ge-Pv and Ge-PPv, respectively. 
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The effect of Fe^^ and Fe^’^ on the Pv to PPv transition in Si-Pv and Ge-Pv is sum¬ 
marized in Fig. [^showing the pressure dependence of enthalpy differences (AH=Hp/),,- 
Hp,.) for the configurations indicated in the previous paragraph. Static (zero tempera¬ 
ture) Pv to PPv transition pressures in iron-free (x=0) MgGeOs and MgSiOs are 43.2 
GPa and 94 GPa, respectively. The Clapeyron slope for the Pv to PPv transitions 
have been previously estimated to be ~ 7.5 MPa/K for MgSiOa ( Tsuchiya et al. | p004l l 
and 7.8 MPa/K for MgGe 03 ( Tsuchiya and Tsuchiya\ 2007| l. Taking these Clapeyron 
slopes approximately into account, transition pressures at ~2500 K should be ~ 62.7 
GPa and ~ 112.8 GPa for MgGe 03 and MgSi 03 , respectively. These results are in 
good agreement with previous experimental measurements and ab initio LDA calcula¬ 
tions ^Murakami et al.\\2QQA\^irose et ^ |2005{ |75Mc/t/yQ et al.\\2Q0A\^suchiya and] 
\Tsuchiy^\2007\ . 

A summary of transition pressure shifts caused by iron in silicate and germanate 
phases is indicated in (Fig. Table [^. As indicated, Fe^^ decreases the Pv to PPv 
transition pressure in both MgGe 03 and MgSi 03 . These observations are consistent 
with previous experimental ( |Mao ef oT 2004[ Shieh et al. 2006| l and ab initio studies 
( Caracas et al. 2005 1 Stackhouse et al.\ 2006| . Flowever, in the case Fe^^, the Pv to 
PPv transition pressure increases noticeably. To the best of our knowledge, this fact 
has not been verified experimentally yet. 

To understand these opposite effects of Fe^'^ and Fe^^ substitution on the Pv to PPv 
transition, the volume difference between Si-Pv and Si-PPv phase (AV=Vpp,, -Vp„) is 
shown in Fig. |^a). With respect to iron-free case, the volume difference is smaller and 
larger in magnitude with Fe^^ and Fe^’^ substitution, respectively. These smaller and 
larger volume difference will contribute to enthalpy difference between Pv and PPv 
phase in such a way that will decrease the trasniton pressure with Fe^^ substitution 
while it will increase when Fe^’^ is substituted. This can be seen explicitly in Fig. [^b 
and c), where pressure dependence of energy (AE) and PxV term [A(PV)] of relative 
enthalpy (AH) are shown separately. Transition pressure points, where AE and A(PV) 
terms cancel each other, are shown by vertical dashed lines. 
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5. Conclusions 


Using ab initio DFT+U calculations, we have undertaken a comparative investiga¬ 
tion of state changes of ferrous (Fe^^) and ferric (Fe^^) iron in MgSi 03 and MgGeOa 
perovskite (Pv) and post-perovskite (PPv) phases. The goal was to verify the possi¬ 
bility of using MgGe 03 as a low pressure analog of MgSi 03 for better understanding 
of spin state changes in iron. Unfortunately, owing to the larger ionic radius of ger¬ 
manium, Fe-O inter-atomic distances in the germanate phases are larger than in the 
silicate phases, allowing the germanates to sustain higher pressures than analog sili¬ 
cates before state changes in iron take place. It has been clearly shown that such state 
changes, whether the displacement transition (low- to high-QS) of Fe^^ in the A-site 
or the HS to LS state change in Fe^^ in the B-site, take place at approximately the 
same average Fe-O, (Fe-O), distance irrespective of the chemistry or the exchange and 
correlation functional used in the calculations. 

This encompassing investigation also allowed us to address the effect of Fe^^ and 
Fe^^ substitutions on the Pv to PPv transition. Despite the simplicity of these static 
calculations, previously measured and calculated effect of Fe^^ on this transition in 
MgSiOs is well reproduced: Fe^^ substitution decreases the transition pressure, whether 
in germanate or silicate phases. However, we report the opposite effect for Fe^^ substi¬ 
tution in the B-site. This transition should involve the LS-B state in both Pv and PPv 
silicate phases but HS-B and LS-B states in Pv and PPv germanates, respectively. For 
both types of transitions the effect is the same: coupled substitution of Fe^^ increases 
the Pv to PPv transition in both materials. This observation has not been validated by 
experiments yet and should be relevant for better understanding the D“ discontinuity 
and the nature of the deep lower mantle. 
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Table 1: Self-consistent Hubbard U.yc(eV) for fen'ous and ferric iron calculated by\Hsu et fl/.|j2010||2011| in 
Pv and hy \Yu et i37^j2Q12^ in PPv structures. _ 





Fe^+ 



High-QS Low-QS 

A(HS) 

B(HS) 

B(LS) 

Pv 

3.1 3.1 

3.7 

3.3 

4.9 

PPv 

2.9 

4.0 

3.5 

5.6 


Table 2: Spin state transition pressures, Pf(GPa), in iron-bearing MgGe 03 and MgSi 03 calculated within 
LDA-I-U.?c and GGA+U^c methods. 



Low-QS to High-QS in -bearing Pv phase 

B(HS) to B(LS) in Fe^^-bearing Pv and PPv phases 



Pv 

PPv 


LDAh-U sc 

GGA+Ujc 

LDAh-TJ^c GGA-hTJ^c 

LDA-bU sc 

GGA+Ujc 

MgGe 03 

28.5 

40.5 

56.5 85.0 

49.5 

82.5 

MgSi 03 

9.5 

22.5 

41.0 69.5 

28.0 

59.5 


Table 3: Perovskite to post-perovskite transition pressures, Pc(GPa), calculated within LDA+U.yc methods. 
X represents the iron concentration. 

X 0 0.125 (Fe2+) 0.125 (Fe^+) 

MgGeOs 432 405 525 

MgSiOs 94.0 88.5 104.0 
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Figure 1: (Color Online) (a) Atomic structure around low-QS and high-QS Fe^'*'-bearing MgSi 03 in Pv phase 
at ambient pressure. Fe, Mg, O, and Si are represented as purple, orange, red, and blue sphere, respectively. 
MgGeOs-Pv structures are similar and not shown here, (b) Enthalpy difference (AH=H/;/g/j_g 5 -l^hw-QS ) for 
Fe^'^-bearing Si-Pv (blue) and Ge-Pv (red). Lower and upper bounds correspond to LDA+Usc and GGA+U.yc 
results, respectively. 
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Figure 2: (Color Online) Atomic environment around iron in -bearing Si-Pv and Ge-Pv phases, (a) Fe- 
O distances (in A) in Si-Pv, and Ge-Pv structures at ambient pressure are shown next to the oxygen atoms, 
(b) Pressure dependence of average Fe-0 distance. Vertical lines represent the low-QS to high-QS transition. 
Lower and upper bounds correspond to LDA+U. 5 C and GGA+U.yc results, respectively. 
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Figure 3: (Color Online) (a) Atomic structure of Fe^'^-bearing MgSiOs Pv and PPv phase. Structures of Pv 
and PPv phases of Fe^'''-bearing MgGe 03 are similar and not shown here. Fe, Mg, O, and Si are represented 
as purple, orange, red, and blue sphere, respectively. Enthalpy difference (AH=H 5 (/^)-HB(y^ 5 )) for Fe'^'^- 
bearing MgSiOs (blue) and MgGeOs (red) in (b) Pv and (c) PPv phases. Lower and upper bounds correspond 
to LDA+U.VC and GGA+U^c results, respectively. 
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(a) 



(b) 



Figure 4: (Color Online) Atomic environment around Fe^'^ in Si-Pv and Ge-Pv phases, (a) Zero pressure 
Fe-0 bond lengths (A) in Si-Pv and in Ge-Pv obtained with LDA+U^c- (b) Pressure dependence of average 
Fe-0 bond lengths for MgSiOs (blue) and MgGeOs (red) phases. Vertical lines represent the high- (HS) to 
low-spin (LS) transitions in the B-site. Lower and upper bounds correspond to LDA-I-U^c and GGA-I-U.?c 
results, respectively. 


21 














30 40 50 60 70 80 90 

Pressure (GPa) 


Figure 5: (Color Online) Enthalpy difference (AH=H/>pv-H/>v) between Fe^^-bearing Ge-Pv and Ge-PPv 
phase using LDA+Ujc method with different possible spin-state combinations of Fe^+ at B-site. (a) Pv(FlS) 
and PPv(HS), (b) Pv(HS) and PPv(LS), (c) Pv(LS) and PPv(HS), (d) Pv(LS) and PPv(LS). 
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Figure 6: (Color Online) Enthalpy difference (AFt=Fl/>/>,,-Ft/>,,) obtained using LDA+U.,c for (a) iron free, (b) 
Fe^"'^-, and (c) Fe^^-bearing Si-Pv and Ge-Pv. 
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Figure 7: (Color Online) Difference in (a) volume {A\=YppyWpy), (b) internal energy (AE=E/>Pv-Epv). and 
(c) PV term (A(PV)=(PV)/-/>v-(PV)/>v) of the enthalpy for iron-free (red), Fe-+-beai'ing (blue), and Fe^+- 
bearing (green) MgSiOs. Results represent LDA-hU^c calculation. Vertical dashed lines represent Pv to PPv 
transition point. 
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